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Small cyclic variations in the frequencies of acoustic modes are expected to be a common phenomenon in solar-like 
pulsators, as a result of stellar magnetic activity cycles. The frequency variations observed throughout the solar and stellar 
cycles contain information about structural changes that take place inside the stars as well as about variations in magnetic 
field structure and intensity. The task of inferring and disentangling that information is, however, not a trivial one. In the 
sun and solar-like pulsators, the direct effect of the magnetic field on the oscillations might be significantly important 
in regions of strong magnetic field (such as solar- / stellar-spots), where the Lorentz force can be comparable to the 
gas-pressure gradient. Our aim is to determine the sun- / stellar-spots effect on the oscillation frequencies and attempt 
to understand if this effect contributes strongly to the frequency changes observed along the magnetic cycle. The total 
contribution of the spots to the frequency shifts results from a combination of direct and indirect effects of the magnetic 
field on the oscillations. In this first work we considered only the indirect effect associated with changes in the stratification 
within the starspot. Based on the solution of the wave equation and the variational principle we estimated the impact of 
these stratification changes on the oscillation frequencies of global modes in the sun and found that the induced frequency 
shifts are about two orders of magnitude smaller than the frequency shifts observed over the solar cycle. 
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1 Introduction 

The level of magnetic activity in stars varies periodically 
with time. The best example of this is provided by the sun 
and the most obvious sign of the solar activity is the peri- 
odic variation of the number of sunspots over the so-called 
solar cycle. The magnetic activity in the sun and other stars 
also influences the properties of the global waves that prop- 
agate within them. In the sun, the variations of the oscilla- 
tion properties along the solar cycle are well documented 
and studied (e.g. Libbrecht & Woodard 1990; Chaplin et al. 
2007): during the solar maximum the frequencies of acous- 
tic modes increase and the amplitudes decrease. Frequency 
shifts were also discovered in a CoRoT star (Garcia et al. 
2010) and are expected to be found in stars observed by Ke- 
pler. 

Sunspots induce direct and indirect effects on the oscil- 
lation properties. The total contribution combines these two 
components. The direct effect results from the action of the 
Lorentz force on the perturbations and the indirect effect is 
due to the magnetically induced changes in the stratifica- 
tion within the spot and due to the effect of the magnetic 
field on the convective flows, which affects the efficiency of 
convection and the turbulent velocities. 

In this paper we describe a possible approach to calcu- 
late some of these effects and present the results for the indi- 
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rect contribution associated with changes to the star's strat- 
ification within starspots. In the following section (Section 
|2l, we present the wave equation and the method adopted 
to solve it. To solve the wave equation we need a starspot 
model. The model used in this work is presented in Section 
[3] Finally, we show our results for the case of the sun in 
Section|ll 

2 Wave Equation 

In this work we consider only the indirect effect of the mag- 
netic field on the oscillations, in particular that associated 
with the stratification changes within the starspot, as com- 
pared with the spherically symmetric case. 

Applying the linear perturbation theory and neglecting 
the perturbation to the gravitational potential and the direct 
effect of the magnetic field on the oscillations (i.e., neglect- 
ing the terms coming from the Lorentz force), we can de- 
rive, in the adiabatic case, the following wave equation 

where £ is the displacement around the equilibrium posi- 
tion, t the time, 7 the first adiabatic exponent, p the density, 
p the pressure, $ the gravitational potential and the index 
'o' denotes the equilibrium quantities. The indirect effect 
of the magnetic field is seen in this equation in the terms 
involving p and p . These quantities differ from the form 
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they take in the spherically symmetric case, depending on 
depth and on the distance from the spot's axis according to 
the model presented in Section[3] 

To solve Equation ([T} we used cylindrical coordinates 
(r, (f>, z) and took advantage of the axial symmetry of the 
spot model (Section 0). In practice we followed the proce- 
dure described below: 



On the outer boundary the combined model merges into 
a quiet sun nonmagnetic model (Model S of Christensen- 
Dalsgaard et al. (1996)) matched to the VAL-C model of the 
solar chromosphere (Vernazza et al. 1981). Since we want 
to integrate the wave equation to layers that are deeper than 
the spot's base, we join the pressure and density profiles of 
the spot into those of the Model S at the lower boundary. 



1. We started by neglecting the horizontal derivatives of 
the displacement by assuming that the displacement 
varies faster in the vertical direction. We also neglected 
the horizontal derivatives of pressure and density, solv- 
ing the equation locally at each distance from the spot's 
axis (r). In this way our problem was reduced to a sec- 
ond order total differential equation for the vertical dis- 
placement. This equation was solved in the region of 
the spot, down to a depth of - 40 Mm (thus, beyond the 
spot's base), by applying the boundary condition that the 
Lagrangian pressure perturbation vanishes at the outer- 
most point of our model, which implies (V • £) = 0, and 
a normalization condition at the surface. 

2. Since the pressure and the density stratifications within 
the spot depend on r, the solution found in point 1. is 
also dependent on r. Thus, we took the approximate so- 
lution derived in 1 . and computed the horizontal deriva- 
tives of the displacement. With this at hand we consid- 
ered again Equation ([T}, now taking into account also 
the horizontal derivatives of pressure and density, and 
derived a new solution for the vertical component of the 
displacement, as well as the solution for the horizontal 
component. 

3. Finally, with the vertical and horizontal components at 
hand from the previous iteration we calculated their hor- 
izontal derivatives and fed them back into the original 
equation for a new solution. This last procedure was re- 
peated until the horizontal and vertical components of 
the displacement remained unchanged. 

In practice, the final solution for the vertical component 
of the displacement differed only very little from the first 
solution, thus confirming that the terms involving the hori- 
zontal derivatives have a relatively small impact on the so- 
lutions for the vertical displacement component. 



3 Starspot Model 

In order to determine the contribution to the frequency 
shifts of the regions of intense magnetic field, we need a 
starspot model in magnetohydrostatic equilibrium. We use 
the sunspot model of Khomenko & Collados (2008) that is 
constructed in cylindrical coordinates, with B, p D and p 
dependent on the distance from the spot's axis, r, and on 
the height measured from the field-free photosphere, z. This 
model combines two other models: the approach of Pizzo 
(1986) for the photosphere and the model of Low (1980) 
that is more appropriate in deep layers. 




Fig.l 

model. 



Magnetic field lines of Khomenko & Collados 



The magnetic field topology of the sunspot model is 
shown in Figure[T[ where z = —10 Mm is the sunspot base 
and z — is the photosphere (z and r are the cylindrical 
coordinates). 

We can see through the pressure and density profiles 
(Figures [2] and [3]l that the Khomenko & Collados model 
does not merge perfectly to the Model S at large distances 
from the spot's axis. Since we are interested in investigat- 
ing the impact of the sunspot structure on the oscillations 
we will thus take as the reference magnetic -free model the 
model of Khomenko and Collados at the spot's outer bound- 
ary (r = 40 Mm). Considering Model S as the reference 
would introduce an artificial difference between the spot 
and non-spot solutions, which would be associated with the 




Fig. 2 Pressure profile of Model S and of Khomenko & 
Collados model for three values of the distance r from the 
sunspot's axis: at the spot's outer boundary (40 Mm), umbra 
(10 Mm) and at the spot's axis (0 Mm). The larger square 
shows a zoom of the region inside the smaller square. 
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Fig. 3 Density profile of Model S and of Khomenko & 
Collados model for three values of the distance r from the 
sunspot's axis: at the spot's outer boundary (40 Mm), umbra 
(10 Mm) and at the spot's axis (0 Mm). The larger square 
shows a zoom of the region inside the smaller square. 

imperfect matching between the stratification at the spot's 
edge and the stratification of the quiet sun. 

4 Results of the Indirect Effect of Sunspots 

In order to evaluate the impact on the oscillations of the 
structural variation within the spot, we computed the solu- 
tions of Equation (fTh for a mode of I = and several fre- 
quencies. In FigureBlwe show the solutions for the radial 
displacement, found following the procedure described in 
Section [2j at different distances from the sunspot's axis for 
v = 3711 /iHz. We emphasize that the differences in these 
solutions result only from the differences in the stratifica- 
tion. No direct effect of the magnetic field on the oscilla- 
tions is considered here. 



o 

cos(e) 

Fig. 5 Phase difference between the reference solution (at 
40 Mm) and the solutions at different distances to the spot's 
axis as a function of the cosine of the photospheric magnetic 
field inclination angle. 

where z* is the depth at the spot's base, 6 the phase and k is 
the acoustic wavenumber 

2 2 

, 2 = ^, O) 

where uj is the angular frequency, uj c is the critical acous- 
tic frequency and c is the adiabatic sound speed. The phase 
difference between the reference solution at 40 Mm and the 
solutions for different distances from the spot's axis is illus- 
trated in Figure [5] as a function of the cosine of 6, where 9 
is the inclination angle of the magnetic field at the photo- 
sphere relatively to the vertical axis. As expected, the phase 
difference changes smoothly from zero, as we approach the 
spot's axis, becoming significant only in the innermost re- 
gions where the pressure and density stratifications are sig- 
nificantly different from their quiet sun counterparts. 



; r=0Mm 
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Fig. 4 Scaled displacement solution (v = 3711 fiHz, I = 
0) for the non-magnetic case at the spot's outer boundary 
(40 Mm), umbra (10 Mm) and at the spot's axis (0 Mm). 



With the vertical component of the solution at hand, we 
compute the phase by matching it to the asymptotic solution 
(Gough 1993) at each r 



cos 



ndz + S(z* , r) 



(2) 



4.1 Induced Frequency Shifts 

The phase difference computed in Section [4] can be used to 
compute the frequency shifts of global modes by following 
the approach of Cunha and Gough (2000). Using a varia- 
tional principle, these authors have shown that the angular 
frequency shift is related to the phase difference by the ex- 
pression 

A - ^ (4) 



Wo ujl r? c- 2 

where AS is the integral phase shift 



A5 



(5) 



In Equation Q R* is the distance from the stellar center to 
the base of the starspot, (r\ 9, </>) are the spherical coordi- 
nates, Y" 1 is a spherical harmonic of degree I and order m, 
r l is the lower turning point of the mode and the index V 
indicates the equilibrium quantities. 
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The modes that we are considering have I = and m 
I i. hence, Y Q ° = 1/ (2y / Tr) and the normalization factor is 

/ / {Yl n fsm6d6d(t> = l. 
Jo Jo 

The phase difference A<5 is different from zero only in- 
side the sunspots. The total effect on the frequencies can 
thus be computed as the superposition of the contributions 
from each sunspot. Since it is impossible to model each 
sunspot separately we will instead adopt the AS computed 
for the sunspot considered in this work as a proxy, assuming 
all sunspots to be alike. 

For a given sunspot, if we centre the coordinate system 
on the spot's axis, the phase difference is independent of (f> 
and we find 

Aw ^ 1 f ma * A<5 (0) sin BdB 
uj ~ 2w2 ff c -2 K -i dr i 

where 9 max coincides with the edge of the spot. 

The integral in the denominator extends only from the 
centre to the base of the sunspot. In that region spherical 
symmetry is maintained and the integral is best computed 
in spherical coordinates. To compute the integral in the nu- 
merator, we convert to cylindrical coordinates. Noting that 

the ratio between r and R* is equal to the sine of 6, which 

r 

is a very small angle, we write 6 ~ — = f and find, for a 

R* 

given sunspot, 

where f max is the ratio between r at the spot's outer bound- 
ary of the penumbra and R* . 

Assuming that all spots are alike, and keeping in mind 
that we are considering modes of degree I = 0, for which 
the weighing by the spherical harmonic is independent of 
the spot's position, we find that the contribution to the fre- 
quency shifts of a number A*" of sunspots combined is 

Ai^ N Jo ma *^rdf 

For the case of sun, we estimated N as the ratio between 
the daily sunspots' area (yearly averaged) at solar maximum 
and the area covered by our spot. Since the spot used as 
a proxy in our work is large (40 Mm from the axis to the 
edge), N estimated in this way is only 3. 

For v = 3711fiHz we found, using Equation (8), a fre- 
quency shift of 2 x 1CP 3 fiHz. In Figure [6] we show the 
frequency shifts computed for a range of frequencies. 
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Fig. 6 Frequency shifts associated with the indirect effect 
of sunspots considered in this work, as a function of fre- 
quency. 

5 Conclusions 

We have computed the frequency shifts that result from 
comparing a spherically symmetric model of the sun and 
an otherwise similar model but with sunspots. We did not 
take into account the direct effect of the magnetic field on 
the oscillations, nor its effect on the local dynamics. The 
only effect considered in the present work is that resulting 
from the changes in the stratification within the spots. 

The values of the frequency shifts found are almost two 
orders of magnitude smaller than the observed frequency 
shifts between solar maximum and minimum. This indi- 
cates that the stratification changes within sunspots do not, 
by themselves, contribute significantly to the overall fre- 
quency shifts along the solar cycle. Since the main con- 
tribution of this indirect effect comes from the umbra, the 
result might increase slightly if we consider the more realis- 
tic case of a greater number of smaller spots. Nevertheless, 
from the present calculation we may already expect this in- 
direct magnetic field contribution to be always significantly 
smaller that the observed frequency shifts. 

Our future objective is to obtain the magnetic solutions 
and compute the induced frequency shifts in order to esti- 
mate the direct effect of starspots to the observed frequen- 
cies. From the experience gained in this work we expect 
that the contribution associated with the direct effect of the 
magnetic field will be comparable to the observed frequency 
shifts only if the phase shift associated with the direct effect 
is of order 1 and extends beyond the umbra. 

Acknowledgements. The authors wish to thank FCT/MCTES, 
Portugal, for financial support through the project PTDC/CTE- 
AST/98754/2008. MC is partially funded by POPH/FSE (EC). We 
are also very grateful to European Science Foundation (ESF) for 
the travel support and the grant for the conference fee. 

References 

Chaplin, W.J., Elsworth, Y., Miller, B.A., Verner, G.A.: 2007, The 
Astrophysical Journal 659, 1749 

© 2012 WILEY-VCH Verlag GmbH & Co. KGaA, Wdnheim 



1036 



A.R.G. Santos, M.S. Cunha & J.J.G. Lima: Asteroseismology and Magnetic Cycles 



Christensen-Dalsgaard, J., Dappen, W., Ajukov, S.V., Anderson, 
E.R., Antia, H.M., Basu, S., Baturin, V.A. Berthomieu, G., 
Chaboyer, B., Chitre, S.M., et al.: 1996, Science 272, 1286 

Cunha, M.S., Gough, D.: 2000, Monthly Notices of the Royal As- 
tronomical Society 319, 1020 

Garcia, R.A., Mathur, S., Ballot, J., Regulo, C, Metcalfe, T.S., 
Baglin, A.: 2010, Science 329, 1032 

Gough, D.O.: 1993, Zahn, K.-R, Zinn-Justin, J., Amesterdam, 399 

Khomenko, E., Collados, M.: 2008, The Astrophysical Journal 
689, 1379 

Libbrecht, K.G., Woodard, M.F.: 1990, Nature 345, 779 
Low, B.C.: 1980, Solar Physics 67, 57 
Pizzo, V.J.: 1986, The Astrophycal Jornal 302, 785 
Vernazza, J.E., Kosovichevm, A.G: 2003, The Astrophycal Jornal 
591,446 



© 2012 WILEY- VCH Verlag GmbH & Co. KGaA, Weinheim 



www.an-journaI.org 



